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Spin Polarization of the Low Density 3D Electron Gas 

F. H. Zong, C. Lin, and D. M. Ceperley 
Dept. of Physics and NCSA, University of Illinois at Urbana- Champaign, Urbana, IL 61801 

To determine the state of spin polarization of the 3D electron gas at very low densities and 
zero temperature, we calculate the energy versus spin polarization using Diffusion Quantum Monte 
Carlo methods with backflow wavefunctions and twist averaged boundary conditions. We find a 
second order phase transition to a partially polarized phase at r s ~ 50 ± 2. The magnetic transition 
temperature is estimated using an effective mean field method, the Stoner model. 

PACS Numbers:71.10.Ca, 71.10.Hf, 75.10.Lp, 05.50.Fk 



i 

CO 



I 

C 

o 
o 



> 

On 
CO 
CO 

in 
o 

CN 
O 

03 



O 

o 



X 



I. INTRODUCTION 

The three dimensional homogeneous electron gas, also 
known as the fermion one component plasma or jellium, 
is one of the simplest realistic models in which electron 
correlation plays an important role. Despite years of ac- 
tive research, the properties of thermodynamic phases 
of the electron gas are still not known at intermediate 
densities.EI In this paper, we study the spin polarization 
phase transition of the three dimensional electron gas at 
zero temperature with recently improved quantum Monte 
Carlo methods. 

There has been recent interest in the low density phases 
spurred by the observation of a ferromagnetic siate in 
calcium hexaboride (CaB@) doped with lanthiumo. The 
magnetic moment corresponds to roughly 10% of the 
doping density. The temperatures (600K) and densities 
(7 x 10 19 /cm 3 ) of this transition are in rough agreement 
with theppredicted transition in the homogeneous elec- 
tron gasJa. However, to make a detailed comparison, it 
is necessary to correct for band effects. For example, 
conduction electrons are located at the X-point of the 
cubic band structure and thus have a six-fold degener- 
acy. The effective mass of electrons at this point and 
the dielectric constant are also changed significantly from 
their vacuum values B These effects cast doubt on the vi- 
ability of the electron gas model to explain the observed 
phenomena. Excitonic models have been proposed to ex- 
plain the ferromagnetismEX Whatever the interpretation 
of ferromagnetism in CclBq, the determination of the po- 
larization energy of the electron gas is important problem 
because of the importance of the model. 

The ground state properties of the electron gas are en- 
tirely determined by the density parameter r s = a/ag 
where Airpa 3 /3 = 1 and ao is the bohr radius, possibly 
changed from its vacuum value by band effects. In effec- 
tive Rydbergs, the Hamiltonian is : 



1 N 9 

j,2 ' 4 1 y° L. 4 

s i=l s i<j 



1 



const. 



(1) 



Note that the kinetic energy scales as \/r 2 s and the po- 
tential energy scales as l/r s so that for small r s (high 
electronic density), the kinetic energy dominates, and the 
electrons behave like ideal gas; in the limit of large r s , 



the potential energy dominates and the electrons crystal- 
lize into a_Wigner crystaLB There is a first order freezing 
transition^ at r s w 100. 

Considering now the spin degrees of freedom, at small 
r s , electrons fill the Fermi sea with equal number of 
up spin and down spin electrons to minimize the to- 
tal kinetic energy and thus the total energy; the sys- 
tem is in the paramagnetic state. As the density de- 
creases and before the freezing transition, there is a pos- 
sibility that the electrons become partially or totally 
polarized (ferromagnetic). The spin polarization is de- 
fined as C = Wi -Ni\ /N, where N-\ and N t are the 
number of up and down spin electrons, respectively and 
N = Ni + iVj,. For paramagnetic phase £ = and for 
ferromagnetic phase £ = 1. 

This polarization transition was suggested by BlochQ 
who studied the polarized electronic state within the 
Hartree-Fock (HF) approximation. He found the fer- 
romagnetic state favored over paramagnetic state for 
r s > 5.45, almost within the density of electrons in met- 
als. However, HF is not accurate for r s > 0. 

More accurate energies became available with the de- 
velopment of Monte Carlo methods for many- fermion sys- 
tems. Ceperleytj using Variational Monte Carlo with a 
Slater- Jastrow trial function determined that the tran- 
sition between the polarized and unpolarized phase oc- 
cured at r s = 26 ± 5. Using a^more accurate method, 
diffusion Monte Carlo (DMC),H it was estimated that 
the polarized fluid phase is stable at r s = 75 ± 5. An 
extension to this workEI found the £ = 0.5 partially po- 
larized fluid becomes stable at roughly r s w 20 and the 
completely polarized state never stable. 

Recently Ortiz et al.tl applied similar methodsETJ to 
much larger systems (N < 1930) in order to reduce 
the finite-size error. They concluded the transition from 
the paramagnetic to ferromagnetic transition is a con- 
tinuous transition, occurring over the density range of 
20 ± 5 < r s < 40 ± 5, with a fully polarized state at 
r s > 40. 

Due to the very small energy differences between states 
with different polarizations, systematic errors greatly af- 
fect the QMC results. Recent progress in the quantum 
simulation methods, makes it possible to reduce these er- 
rors. Kwon et aLEa found that a wavefunction incorpo- 
rating back-flow and three-body (BF-3B) terms provides 



a more accurate description: they obtained a significantly 
lower variational and fixed-node energy. In another ad-, 
vance of technique, twist-averaged boundary conditionsO 
(TA) have been shown to reduce the finite-size error by 
more than an order of magnitude, allowing one to obtain 
results close to the thermodynamic limit using results for 
small values of N. In this paper, we apply these improved 
methods to the polarization transition in the three di- 
mensional electron gas. We first describe the simulation 
method, and then, the results. 



II. METHODS 

The most accurate Quantum Monte Carlo (QMC) 
methodE-l at zero temperature is Projector or Diffusion 
Monte Carlo (DMC): one starts with a trial function and 
uses exp(— tH) to project out the ground state using 
a branching random walk. Fermi statistics pose a sig- 
nificant problem for the projection method, since exact 
fcrmion methods such as transient estimate or release- 
node QMC suffer an exponential loss of efficiency for 
large numbers of particles. For this reason, the fixed- 
node approximation is normally used, obtaining the best 
upper bound to the energy consistent with an assumed 
sign of the wave function. The generalization of the fixed- 
node method to treat complex-valued trial functions is 
known as the fixed-phase approximation.!!-! 

In the simpler, but less accurate Variational Monte 
Carlo method (VMC), one assumes an analytic form 
for a trial function ^>t{R) and samples l 1 !^.??)! 2 us- 
ing a random walk. An upper bound to the exact 
ground state energy is the average of the local energy 
E L {R) = ^ T (R)- 1 H^ T {R) over the random walk. 

The trial wavefunction plays a very important role in 
these two methods. With a better trial wavefunction, 
not only is the variational energy lower and closer to the 
exact energy, but also the variance of the local energy 
is smaller so that it takes less computer time to reach 
desired accuracy level. The trial wavefunction is also 
important in fixed-phase DMC because the solution is 
assumed to have the same phase as the trial function. 
One then solves for the modulus. This implies that the 
DMC energy lies above the exact ground state energy 
by an amount proportional to the mean squared differ- 
ence of the phase of the trial function from the exact 
phase. As this is the only uncontrolled approximation, 
it is important to carefully optimize the assumed trial 
wavefunction. 

For a homogeneous system, the non-interacting (NI) 
wavefunction consists of a Slater determinant of single- 
electron plane waves orbitals. To incorporate electron 
correlation, one multiplies the NI wavefunction by a pair 
wavefunction, obtaining the so-called Slater- Jastrow (SJ) 
form. To construct a better trial function, one incor- 
porates backflow and three-body effects. 13 The particle 
coordinates appearing in the determinant become quasi- 



particle coordinates: 

N 

Xi = r i + ^7 ? (r ij )(r i -r j ). (2) 

The Slater determinant is then D = det(e' km ' Xn ) where 
rj(rij) is a function to be optimized. Then the backflow- 
three body wavefunction (BF-3B) is: 
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where 

N 

Gi = £ffoOfo-ri) ( 4 ) 

and 

u(r) = u RPA (r) - AtC 2 (r)r 2 + 7 (r). (5) 

Here and are determinants for the up spin and 
down spin electrons, £(r) is the three-body correlation 
function, and u is the Jastrow correlation fuaction. For 
the electron gas an accurate analytic form,__l u RPA (r), 
has as low an energyQ as those with optimized parame- 
ters. In the presence of three-body correlation, the RPA 
two-body term is supplemented with an extra Gaussian 
function j(r). Please refer to Kwon et al.t-l for further 
details concerning this wavefunction. We used optimized 
Ewald sumst-l both for the potential and for the corre- 
lation factor so as to have the correct long wavelength 
behavior. 

Though the computational cost for BF-3B wavefunc- 
tion is somewhat greater than the simple SJ wavefunc- 
tion and there is the added cost of optimization for BF- 
3B wavefunction, we found accurate trial wavefunctions 
crucial to compute the small energy differences between 
different polarization states. We optimized the parame- 
ters by minimizing a combination of the energy and the 
variance for each density and polarization. Figure (1) 
shows the energy vs. polarization at r s = 50 using dif- 
ferent trial functions and simulation methods. The SJ 
trial function with VMC has the highest energy for all 
polarizations and at this level of accuracy finds the fully 
polarized phase to_.be stable, in agreement with earlier 
VMC calculations J_ However, using the best BF-3B trial 
function, the variational energies are lowered significantly 
with the unpolarized energy dropping more than the po- 
larized case so that the polarized phase is no longer sta- 
ble. DMC calculations confirm this result. Note that the 
DMC energies determined using the NI phases (or nodes) 
give energies lower than the BF-3B variational energies 
confirming the importance of accurate DMC calculations. 
The use of BF-3B wavefunctions with DMC leads to the 
lowest ground state energies, hopefully, very close to the 
exact energy. 
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FIG. 1. Energy versus spin polarization at r s — 50 for 54 
electrons using TA with 10 3 twist values. Compared are calcu- 
lations with SJ and BF-3B wavefunctions and with two QMC 
methods: VMC and DMC. 

After the effect of the nodes, the dependence of the 
energy on the number of electrons is the largest system- 
atic error. Within periodic boundary conditions (PBC), 
the phase picked up by the wavefunction as a particle 
makes a circuit across the unit cell, is arbitrary. General 
boundary conditions are: 



*(ri+L,r 2 ,...) = e 4e *(r 1 ,r 2 ,. 



(6) 



where L is a lattice vector of the supercell. If the twist an- 
gle 9 is averaged over, most single-particle finite-size ef- 
fects arising from shell effects in filling the plane wave or- 
bitals, are eliminated. This is particularly advantageous 
for polarization calculations since shell effects dominate 
the polarization energy. The extra effort in integrating 
over the twist angles is minimal, since the various calcu- 
lations all serve to reduce the final variance of the com- 
puted properties. The effect of-iboundary conditions is 
examined in detail in Lin et aZJla 

There is a further size effect in the calculation of the 
potential energy due to a charge interacting with its cor- 
relation hole in peighboring supercells as shown in Fig. 
(6) of Lin et alt3. To correct this, we fit the energies 
versus N using the expansion: 
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E N - E x + — + — + 



(7) 



For unpolarized systems, the fitted Eqc agrees with the 
previous non-twist averaged (PBC) result determined us- 
ing extrapolations based on fermi liquid theory (FLT). As 
shown in Fig. g, the correlation hole is only weakly de- 
pendent on spin polarization at low density: the peak of 
g(r) only changes from 1.175 to 1.190 as the system goes 
from unpolarized to polarized. Hence, the potential size- 
effect hardly changes the spin polarization energy. Fig. 



H shows the polarization energies for N = 54, 108, 162. 
With TA boundary conditions, there is a remarkable in- 
sensitivity to the number of electrons. Even though the 
system size is increased threefold, the change in the en- 
ergy versus polarization is almost undetectable. Consid- 
ering only the leading 1/N correction, we estimated -Boo 
with N = 54 and N = 108. 
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FIG. 2. The pair correlation function for several polariza- 
tions at r s — 50 using DMC. The various curves are for 
f = 0, 0.33, 0.67, 1 with the structure increasing with spin 
polarization. 

Also shown in Fig. ||is the [estimate of the polarization 
energy from Ceperley-Alder.tJ In that work, size effects 
were estimated with Fermi liquid corrections. Rather 
than BF-3B wavefunctions, corrections using the poten- 
tially exact, release-node method were used. The re- 
sults with PBC and FN-DMC are in agreement with the 
present FP-DMC calculations. However, the present re- 
sults have an error bar more than an order of magnitude 
smaller than those of Ceperley-Alder, primarily due to 
increased computer performance. 

Though PBC with FLT corrections are adequate for 
unpolarized and fully polarized systems, the precision is 
limited for intermediate polarizations. To estimate finite 
size effects within FLT, one must perform accurate DMC 
simulations for widely varying system sizes. In Ortiz et 
al.p the simulation size varied from 725 < N < 1450. 
Within DMC it is very time-consuming to ensure uni- 
form accuracy independent of particle number, so that 
one typically determines size effects within VMC, using 
the more approximate SJ trial functions. As we have seen 



in Fig. (1), the SJ trial functions are unreliable at low 
densities. TA boundary conditions allow a much better 
way to estimate energies in the thermodynamic limit of 
partially polarized fermi liquids since the number of elec- 
trons can be held fixed as the spin polarization varies. 
Small system sizes, allowing use of more accurate but ex- 
pensive trial functions and even exact fermion methods, 
give precise estimates of the spin polarization energy in 
the thermodynamic limit. 
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FIG. 3. The polarization energy for various sized systems 
at r a — 50 using TA and DMC. (circle 54, square 108, and 
triangle 162). The peint at £ = 1 with the large error bar is 
from Ceperley- Alder Other errors are less than 10 -6 . 
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FIG. 4. Comparison of this work with that of Ortiz et al.0 
The filled diamonds are DMC simulations with TA and BF-3B 
wavefunctions (this work). rjThe filled circles are DMC with 
PBC and SJ wavefunction.EJ All energies are extrpolated to 
the thermodynamic limit. Errors are given in Table I and are 
smaller than the size of the points. 

In Fig. H we compare the total energy for r s = 40 
calculated with DMC and TA and extrapolated to the 
thermodynamic limit with the calculation of much larger 
systems (N = 725) of Ortiz et alB The BF-3B ener- 
gies are lower energy and show a different polarization 
energy: Ortiz's calculation finds that the polarized or 
partially polarized phase is stable at this density, while 
we find the unpolarized phase is stable. This difference 
is due to the backflow correlations in the trial wavefunc- 
tion. Although backflow energies are small, they favor 
the unpolarized state and hence are crucial for accurate 
determination of the polarization transition. 



III. RESULTS 



We carried out computations of the spin polarization 
energies at electronic densities 40 < r s < 100. At each 
density, we performed DMC calculations with N = 54 
and N = 108 electrons using 10 3 twist angles. The time 
step was adjusted so that the DMC acceptance ratio was 
in the range 98% — 99%. Note that when calculating the 
polarization energy, most time step errors will cancel out 
of the polarization energy. Thus systematic errors in the 
polarization energy are much smaller than in the total 
energy. We then extrapolated the energy to the thermo- 
dynamic limit using Eq. (]?]). The energies are given in 
Table I. 
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FIG. 5. The spin polarization energy of the 3DEG times 
r^ 2 in Ry/electron at various densities using a polynomial fit 
to the data in Table I. The density, r s , is denoted on the right 



FIG. 6. The square of the spin polarization versus r s . The 
curves were obtained using fits in the Fig. (4). The line is a 
fit though the points. The value at r s — 40 was obtained by 
extrapolation from physical values of £. 



We then fit the energy versus polarization to a 
quadratic polynomial in £ 2 . The results are shown in fig- 
ure H A polarization transition is evident. At r s = 40, 
the system is still paramagnetic, with the unpolarizcd 
phase stable. As the density decreases, at r s m 50, the 
system becomes unstable with respect to spin fluctua- 
tions. The partially polarized states become stable at 
r s > 60. As the electronic density continues to decrease, 
the fully polarized state has a lower energy with respect 
to unpolarized state at r s > 80, however, we find that 
the partially polarized state has an even lower energy. 



The quoted error bar on the critical density estimates 
the statistical errors, not the systematic errors arising 
from the fixed-phase approximation. The experimental 
and theoretical results on polarized helium give caution 
on placing too much confidence in the estimate of the po- 
larization transition. Even using the accurate optimized 
BF-3B wavefunctions, the magnetic susceptibility in liq- 
uid 3 He does not agree with experiment at low pressure 
and the polarized phase is nearly degenerate with the 
unpolarized phase at the freezing densityEH The present 
results also do not preclude the existence of phases with 
other order parameters such as superfluids, as occurs in 
ground state of liquid 3 He. In fact, it is rather likely 
that the ground state of the electron gas will have such 
a phase at the lowest fluid densities. 



In Fig. |6| is shown the predicted square of the optimal 
polarization versus density. We find that the equilibrium 
polarization is described by £ 2 = (r s — r*)/62 with the 
critical density r* = 50 ±2. As the density decreases, the 
stable state becomes more and more polarized approach 
fully polarized at freezing, r s w 100. Quantum critical 
fluctuations, not present in systems with N < 162, could 
modify the behavior of the spin polarization energy near 
the critical density. 



However, examination of the variance of the trial func- 
tion suggests that the result for the electron gas may be 
more reliable than for liquid 3 He. Shown in fig. ^ is 
the variance of the trial function at r s — 50 for both 
the SF and BF-3B functions. Although the variance of 
the SJ trial function depends on spin polarization, that 
of the BF-3B does not. Such is not the caseO with lia-. 
uid 3 He. Arguments based on variance extrapolationta 
suggest that the DMC calculations with BF-3B phases 
should be more reliable than in liquid 3 He. 
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FIG. 7. The variance of the SJ and BF-3B trial function as 
a function of spin polarization at r s — 50 using TA and VMC 
with N = 54. 



IV. THE PHASE DIAGRAM 

One can use the calculated energies to estimate t 
finite temperature behavior within the Stoner mndel,c3 
arising in the theory of itinerant magnetism.E3 The 
Stoner model differs from Hartree-Fock by replacing the 
Coulomb interaction by a zero range one, a repulsive 
delta function potential: Yli<j 9^( r ij)- One can view 
this approach as the first step to a full Fermi-liquid de- 
scription of the quasiparticle interactions, and use the 
QMC data to determine the strength of those interac- 
tions. One expects that the Slater- Jastrow trial function 
has screened off the long-range interaction leaving only a 
short-ranged spin-dependent term that can be modeled 
by a contact interaction. 

In the Stoner model, the energy is evaluated within 
the mean field (Hartree-Fock) approximation using the 
NI wavefunction. The energy at zero temperature in the 
thermodynamic limit is: 

E oc (1 + C) 5/3 + (1 - C) 5/3 + 0.0545rf(l - C 2 ). (8) 

For gr 2 s < 20.5 the system has an unpolarized ground 
state and for gr 2 s > 24.4 the ground state is ferromag- 
netic. For intermediate couplings, the ground state has a 
partial spin polarization at zero temperature, similar to 
the observed behavior of the electron gas at low density. 

Although the polarizations are qualitatively correct, 
the above functional form does not fit well the DMC data 



{i.e. from Table I.) In addition, assuming g does not have 
a very strong density dependence, the Stoner model pre- 
dicts that the partially polarized density range should be 
quite narrow, from 50 < r s < 54, while as the QMC re- 
sults indicate a much broader density range. Perhaps the 
assumption of a zero-range interaction of quasiparticles 
is too restrictive. 
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FIG. 8. The phase diagram of the electron gas. Conver- 
sion to units of cm and Krwas done using ao = 1.3nm and 
Ry — 250K using estimates^ of the effective mass and the di- 
electric constant of SrB?,. The solid line the mean-field esti- 
mate of the magnetic transition temperature from the Stoner 
model, where the spin interaction is estimated from the zero 
temperature QMC data. The dotted line is the energy differ- 
ence between the unpolarized and partially polarized system. 



Nevertheless, we use this mean field model to make a 
estimate of the transition temperature of the polarized 
phase. Fitting the energy to the function given in Eq. 
(8), we obtain estimates of g and the effective mass of the 
electrons. Then, using the free energy of the ideal spin 
(1/2) partially polarized fermi gas with the additional 
spin interaction, we estimate the temperature at which 
the system becomes polarized by determining when the 

rin-stiffness of the unpolarized system vanishes. Fig-, 
is the estimated phase diagram of the electron gas.EU 
In this diagram, the effective mass and dielectric con- 
stant for &tBq, a closely related material to CaB§, have 
been usedcl to convert to units of K and cm 3 . Note that 
both the temperatures and densities of our calculated 
magnetic transition are four orders of magnitude smaller 
than that found experimentally!] in CaB§. Even assum- 



ing errors because of uncertainties in material properties 
and from the mean field methods to estimate T c , these 
estimates are very difficult to reconcile with experiment, 
apparently ruling out an electron gas model of ferromag- 
netism in this material. Also plotted on the phase dia- 
gram is the energy difference between the partially polar- 
ized fluid and the unpolarized fluid as another estimate 
of the magnetic transition temperature. Finite tempera- 
ture QMC calculations would be desirable to confirm the 
mean field estimate of T c . A rough estimate of the limit 
of stability of the Wign<p,crysta£3 is also shown. 

Tanaka and IchimarucJ have computed the polariza- 
tion phase diagram of the electron gas both at zero and 
non-zero temperature using an integral equation method. 
At zero temperature they obtain a result similar to Or- 
tiz at alu, with a continuous transition at r s w 20 and 
a fully polarized fluid state at a slightly higher density 
r s w 22. Apparently this approach, built on local field 
corrections to free fermion response functions, is biased 
by the initial Hartree-Fock assumption and overempha- 
sizes the tendency for ferromagnetism. 

There is recent work on the low deusity electron gas 
in 2D, studying both the Fermi liquidE-3 and the Wigner 
crystal phase. Using QMC similar techniques in the liq- 
uid phase, a polarized phase is found to be stable be- 
tween 26 < r s < 35, though the energy differences are 
even smaller than in 3D. The partially polarized phase 
is never stable. In the 2D Wigner crystal, path integral 
methodsES were used to derive directly the spin Hamil- 
tonian. It was found that the ground magnetic state 
is a spin liquid though the ferromagnetic state has only 
a slightly higher energy at melting. Analogous calcula- 
tions of the-magnetic phase diagram of the WC in 3D are 
underway.ES 



for the calculation of polarization energies, enabling size- 
converged results with fewer electrons. Using relatively 
small system sizes allows one to use more accurate trial 
wavefunctions and to fully converge the diffusion Monte 
Carlo calculations. We find a second order transition to 
a polarized phase at r s = 50 ± 2. 

In general, methods based on the variational principle 
such as the fixed-node quantum Monte Carlo method for 
many-fermion systems, favor phases with a higher sym- 
metry, in this case the Wigner crystal and polarized fluid, 
over the more complicated unpolarized phase, (effectively 
a two-component mixture of spin up and down electrons.) 
Recall that in HF one has a polarization transition be- 
cause antisymmetry is the only way to correlate electrons, 
however, the correlation is only between like spins; hence, 
there is an instability to polarize the system once the po- 
tential is dominant. But using a SJ wavcfunction both 
like and unlike electrons are highly correlated. Our re- 
sults demonstrate that the SJ wavefunctions still pref- 
erentially favors the "simple" phases, even using DMC 
method. This symmetry argument explains the histori- 
cal tendency for the polarized phase become stable over 
a more and more restricted range of density and tem- 
perature as more accurate methods are used. Our QMC 
results using the BF-3B wavcfunction indicate that there 
is still an instability for spin polarization at very low den- 
sity. Although examination of the variance indicates that 
the polarized and unpolarized BF-3B trial functions are 
equally inaccurate, it is still not clear if our finding of a 
polarization transition is an artifact of the assumed trial 
wavefunction. Calculations with the new methods (TA 
and BF-3B wavefunctions) but with the exact fermion 
methods are desirable to resolve the phase of the elec- 
tron gas at intermediate densities. 



V. CONCLUSION 
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TABLE I. Energy of the 3DEG computed using TA- DMC and extrapolated from TV = 54 and N = 108 with (10 3 twist 
angles). Energies are in Ry/electron. The numbers in parenthesis are standard errors in units of W~ 8 Ry. 
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40 -0.03523748(60) 

50 -0.02889900(62) 

60 -0.02452017(44) 

70 -0.02131429(41) 

75 -0.02001137(35) 

85 -0.01784017(30) 

100 -0.01535357(30) 



-0.02889900(66) 
-0.02451866(51) 
-0.02131381(40) 
-0.02001191(37) 



-0.03523295(67) 
-0.02889962(68) 
-0.02452031(48) 
-0.02131621(39) 
-0.02001376(36) 
-0.01784152(32) 
-0.01535340(30) 
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-0.02451963(50) 
-0.02131716(37) 
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-0.02888835(62) 
-0.02451747(42) 
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-0.01784300(32) 
-0.01535639(26) 



-0.02887542(70) 
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